Ordering temperatures of Ising Spin Glasses 



I. A. Campbell 1 

1 Laboratoire des Verves, Universite Montpellier II, 34095 Montpelher, France 
(Dated: February 2, 2008) 

Exploiting an approach due to Singh and Fisher I show that in the high dimension limit the 
ordering temperature of near neighbour Ising Spin Glasses drops linearly with the kurtosis of the in- 
teraction distribution, in excellent agreement with accurate high temperature series data of Daboul, 
Chang and Aharony. At lower dimensions the linear relation no longer applies strictly but the kur- 
tosis can still be taken to be an appropriate parameter for ranking different systems. I also compare 
the series estimates with simulation and Migdal-Kadanoff estimates where these are available. 
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In order to establish precise values of critical expo- 
nents at a continuous transition, it is essential to first 
obtain accurate and reliable values for the critical tem- 
perature T c . The ordering temperatures T c of Ising spin 
glasses (ISGs) are not expected to be universal but to de- 
pend on the form of the interaction distribution. In ISGs 
with random near neighbour interactions of zero mean 
(Jij) — the interaction distributions are normalized ac- 
cording to the standard convention (J?-) = 1 (where (..) 
represents an average over the interaction distribution 
Jij between pairs of spins). It was suggested on purely 
empirical grounds 1 that T c values should be ranked from 
highest to lowest according to the value of the kurtosis 
of the interaction distribution, 

k = {4)/(4) 2 (i) 

(An alternative definition of kurtosis, sometimes referred 
to as the excess kurtosis, corresponds to K — 3 which 
leads to a zero value for the Gaussian distribution. For 
present purposes it turns out that the the definition in 
Eq. Q is more convenient.) Extensive high temperature 
series calculations by Daboul, Chang and Aharony pro- 
vide accurate values of T c in dimensions d from 8 to 4 
which are consistent with this phenomenological rule in 
each dimension 2 . 

Exploiting an approach due to Singh and Fisher— I 
show rigorously that in the high dimension limit the or- 
dering temperature of near neighbour ISGs drops lin- 
early with the kurtosis of the interaction distribution, 
and find excellent quantitative agreement between the 
Singh-Fisher predictions and the data of Daboul et a? 
in dimensions 7 and 8. At lower dimensions the linear 
relation no longer applies strictly but the trend remains 
and the kurtosis can still be taken to be an appropriate 
parameter for ranking different systems. 

Singh and Fisher— used a 1/d expansion technique 
(whose principle largely pre-dates RGT— ) to evaluate T c 
in ISGs as a function of dimension and of interaction dis- 
tribution. On this approach and for the near neighbour 
ISG with bimodal ± J interactions, they define 



and find that 

K]-°' 5 = (cr) a5 [l - (7/2)/a 2 - (21/2)/a 3 - ...] (3) 

where a = (2d— 1). They summed this series, truncating 
appropriately, to estimate the sequence of T C (±J) values 
for the bimodal ±J interaction in different dimensions. 
They obtained satisfactory agreement with the high tem- 
perature series values which had just been been evaluated 
for this interaction^Sil. More recent high temperature se- 
ries work 2 ^ gives values consistent with the earlier series 
estimates for the bimodal distribution. 

For present purposes the important step of Singh and 
Fisher was the next one : they state that for an arbitrary 
interaction distribution that we will label X, the mean 
over the interaction distribution Jij(X) of 

uj c = (tanh 2 [J y (A)/T c (A)]) (4) 

is the same for all distributions to order 1/er 2 . This is 
an implicit equation for T C (X) once lo c is known. The 
leading deviation was given as 

[lu c (X) - w c (±J)] - w c (±J)[l - 3( M 2 - 1)At 3 )-] (5) 

where 

Mn = (tanh 2n [J y /r c (A)])/(tanh 2 [J y /r c (X)])" (6) 

It turns out that the effect of this deviation term is nu- 
merically rather weak in all dimensions. 

If d is high enough, T c values will also be high ( in d = 8 
T c ~ 3.6 in the examples discussed below) so the whole 
range of [Jy /T C (X)] values can be taken to be small and 
we can evaluate tanh 2 (a;) for small arguments only. Now 

tanh 2 (.x) = x 2 - {2/3)x 4 + (17/45)a; 6 ... (7) 

so keeping only the two leading terms, Eq. becomes 

uj c = (J 2 (A))/T C (A) 2 - (2/3)(4(A))/T c 4 (8) 

As all ISG interaction distributions are normalized such 
that (Jfj(X)) = 1 and in the high d limit uj c is indepen- 
dent of distribution, 



w c (±J) = tanh 2 [J/T c (±J)] 



(2) 



uj c = 1/T C (X) 2 - (2/3)K(X)/(uj c T c (X) 4 ) (9) 
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FIG. 1: (Colour on line) T c estimates for near neighbour inter- 
action ISGs in dimension 8 from high temperature series data, 
Daboul it et al« (black squares), the Singh- Fishes- estimates 
(see text) (green triangles), and the leading approximation 
Eq. IIUL straight line. The points for successive K refer to 
ISGs with bimodal, double triangle, uniform and Gaussian 
distributions from left to right (sea- for the definitions of the 
interaction distributions with the standard normalizations) 



FIG. 2: (Colour on line) T c estimates for near neighbour in- 
teraction ISGs in dimension 7 from high temperature series 
data, Daboul it et al- (black squares), the Singh- Fishes- es- 
timates (see text) (green triangles), and the leading approxi- 
mation Eq. (IIUL straight line. The points refer to ISGs with 
bimodal, double triangle, uniform and Gaussian distributions 
from left to right. 



with K(X) the kurtosis of the distribution X. To the 
same level of approximation the equation can finally be 
rewritten 

T C (X, d) = (lK(d)) a5 [l - K(X)u c (d)/2>] (10) 

This should be exact in the limit of high d and small 
to moderate kurtosis. The critical temperature should 
then vary linearly with the kurtosis of the distribution 
K(X), with a slope which is given by [uj c (d)]°- 5 /3. Lo c {d) 
can be fixed if one T C (X) (generally the bimodal value) 
is known. 

I calibrate the single free parameter Lo c (d) by using the 
Daboul et al value for T C (±J) in each dimension (the 
Singh-Fisher 1/a series T C (±J) values are in fact very 
close to the Daboul et al values), except in dimension 
4 where the error bar on the Daboul et al value is high 
and I calibrate with the best simulation value. The io c 
values used are 0.069877,0.08251,0.1395 and 0.2208 in 
dimensions 8, 7, 5 and 4 respectively. The distributions 
studied by Daboul et al were the bimodal, double trian- 
gle, uniform and Gaussian with K — 1,4/3,9/5 and 3 
respectively. In simulations- the Laplacian (or decreas- 
ing exponential) distribution which has K = 6 was also 
studied. 

In Figures 1 and 2 I show the high temperature series 
T c values^ plotted against kurtosis for dimensions d = 8 
and d = 7. The leading approximation, Eq. (JTUJ, gives 
the straight line which lies very close to the Daboul et al 
data points for the different distributions. Using the full 
Singh-Fisher expressions (Eqs. 10} , JSJl, and ©), eval- 
uated by numerical integration without appealing to the 



small argument approximation gives calculated points 
which are in almost perfect agreement with the Daboul 
et al data in these dimensions. The weak deviation from 
the straight line for K = 3 (the Gaussian distribution) is 
due to the error introduced by truncating Eq. J7J) in the 
wings of the Gaussian. 

By dimension 5, figure 3, the first order approxima- 
tion has become much poorer as could be expected a pri- 
ori because the small argument development is no longer 
valid in this lower dimension, where the T c values are 
lower. There are weak deviations of the Singh-Fisher esti- 
mates (without the small argument approximation) from 
the high temperature series data points, meaning that 
Eq. (JSJ| is no longer accurately obeyed, either because 
higher order terms give non-negligible contributions or 
more fundamentally because the 1 j a expansion becomes 
inaccurate at lower dimensions. As pointed out by Singh 
and Fisher their expansion "is at best asymptotic and 
optimum estimates come by truncating after the small- 
est term"—, a procedure which introduces uncertainties at 
lower dimensions. Nevertheless the general trend of de- 
creasing T c with increasing kurtosis is still clearly present. 

Finally for d = 4 the small argument approximation 
to the Singh-Fisher expressions is certainly inappropriate 
but one can again evaluate the full Singh-Fisher ordering 
temperature estimates through numerical integration. In 
this dimension ordering temperatures have also been esti- 
mated independently by simulations, using various com- 
plementary methods. These include a technique rely- 
ing on consistency between dynamic and static critical 
parameters^, domain wall analysis 10 , and Binder param- 
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FIG. 3: (Colour on line) T c estimates for near neighbour in- 
teraction ISGs in dimension 5 from high temperature series 
data, Daboul it et al- (black squares), the Singh- Fishej- es- 
timates (see text) (green triangles), and the leading approxi- 
mation Eq. 11LIL straight line. The points refer to ISGs with 
bimodal, double triangle, uniform and Gaussian distributions 
from left to right. 



FIG. 4: (Colour on line) T c estimates for near neighbour in- 
teraction ISGs in dimension 4 from high temperature series 
data, Daboul it et al- (black squares), the Singh-Fishei- esti- 
mates (see text) (green triangles), simulation estimates- (red 
circles), and MK estimates with 6* = 2.49- (blue inverted 
triangles). The points refer to ISGs with bimodal, double 
triangle, uniform, Gaussian and Laplacian distributions from 
left to right. 



eter measurements 11,12,13,14 . In the two cases (bimodal 
and Gaussian interaction distributions) where simulation 
results using different methods exist, agreement between 
the various simulation estimates is excellent, see Table 
I. In addition, the T c values have been studied in d = 4 
(and d = 3) using a Migdal-Kadanoff (MK) approach 9 
where the optimum size b* of the MK renormalization 
cell is chosen by intrapolating the number of segments b 
in an MK branch until T C {MK) is equal to T c for the ±J 
distribution. T C {MK) values for the other distributions 
are then estimated with this same b* . For d = 4 the opti- 
mum intrapolated value for b was found to be b* = 2.492 
. The different T c estimates are given in Table 1 and 
are shown in Figure 4 where the same general trend once 
more appears. It can be seen that while the Singh-Fisher 
values and the MK values lie close to but not precisely on 
the series and simulation points, the agreement between 
the simulation estimates and the series values is excellent. 
As the simulation calculations arc necessarily carried out 
on samples of moderate size (a maximum of L = 12 up 
to now in dimension 4 ), a priori there might have been 
a lingering suspicion that the T c estimates from simula- 
tions could somehow be intrinsically in error because of 
corrections to finite size scaling that had not been prop- 
erly taken into account. Because the series estimates are 
not subject to this type of error, the agreement beween 
the two sets of quite independent T c values lays this sus- 
picion to rest. One can have confidence that at least in 
dimension d = 4 the simulation estimates of T C (X) ( and 
in particular those using the static/dynamic consistency 
technique 1 which are the most complete ) are intrinsically 



TABLE I: Comparison of estimates of T c in dimension 4 



distribution 


method 


reference 


T c 


bimodal (A" = 1) 


series 




2.02(3 


bimodal 


series 


2 


1.96(7 


bimodal 


simulation 


14 


2.03(3 


bimodal 


simulation 


10 


2.00(4 


bimodal 


simulation 


11 


2.00(1 


bimodal 


simulation 


1 


1.99(1 


triangles(A = 4/3) 


series 


2 


1.95(1 


uniform(A = 9/5) 


series 


2 


1.90(2 


uniform 


simulation 


1 


1.91(1 


uniform 


MK 


9 


1.87(2 


Gaussian(A = 3) 


series 


2 


1.78(1 


Gaussian 


simulation 


13 


1.80(2 


Gaussian 


simulation 


12 


1.80(2 


Gaussian 


simulation 


1 


1.76(2 


Gaussian 


MK 


!) 


1.75(2 


Laplacian(A = 6) 


simulation 


1 


1.52(1 


Laplacian 


MK 


!) 


1.54(2 



correct and are not polluted by uncontroled corrections 
to finite size scaling. 

Empirically, Table 1 shows that the MK T c values are 
a good indication of the true T c values if the appropriate 
choice of b* has been made by calibrating b using the T c 
for one particular distribution (see the discussion in 9 ). 
In dimension 3 where the series method would require 
a very long set of approximants, and where it turns out 
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that in simulations corrections to finite size scaling are 
much stronger than in dimension 4, it would appear that 
the ratios of T c values for different distributions obtained 
through the MK approach*^ should be rather reliable 
and a very helpful indication for controling simulation 
measurements . 

We can conclude that the key to understanding how T c 
values vary from distribution to distribution in ISGs with 
near neighbour interactions is the Singh-Fisher statement 
that u) c (d) ( Eq. ) is almost distribution independent 
for fixed dimension d. The linear dependence of T C (X) 
on the kurtosis K{X) with a slope [w c (d)]°- 5 /3 becomes 
exact in the high dimension limit and is still a close ap- 
proximation in high but finite dimensions (d = 8 and 7). 
Even at moderate d (d = 5 and 4), the general trend 
still holds so the ranking of distributions by their kurto- 
sis remains physically appropriate in agreement with the 



initial phenomenological suggestion. 1 . ( It is nevertheless 
conceivable that for low dimensions differences between 
distributions could lead to situations such that two dis- 
tributions with identical kurtoses would have marginally 
different T c values. For the moment the number of differ- 
ent distributions that have been studied is very limited 
so this possibility has never been put to the test). In 
dimension d = 4 which is the only dimension where ap- 
propriate data is available for both simulations and series 
measurements, simulation estimates of ordering temper- 
atures are in excellent agreement with the series values 
which confirms that the simulation estimates, in particu- 
lar those based on the consistency between dynamic and 
static critical parameters^, are accurate and reliable. 
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